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We study the ground state phase diagram of ultracold dipolar gases in highly anisotropic traps. 
Starting from a one-dimensional geometry, by ramping down the transverse confinement along one 
direction, the gas reaches various planar distributions of dipoles. At large linear densities, when 
the dipolar gas exhibits a crystal-like phase, critical values of the transverse frequency exist below 
which the configuration exhibits novel transverse patterns. These critical values are found by means 
of a classical theory, and are in full agreement with classical Monte Carlo simulations. The study 
of the quantum system is performed numerically with Monte Carlo techniques and shows that the 
quantum fluctuations smoothen the transition and make it completely disappear in a gas phase. 
These predictions could be experimentally tested and would allow one to reveal the effect of zero- 
point motion on self-organized mesoscopic structures of matter waves, such as the transverse pattern 
of the zigzag chain. 



I. INTRODUCTION 

Ultracold atoms and molecules are an attractive play- 
ground for studying fundamental properties of matter. 
Theoretical and experimental investigations pursue a 
fundamental understanding of the quantum dynamics of 
matter at ultralow temperatures and explore applications 
for quantum technologies, such as quantum metrology 
and information processing p], Q . One relevant issue in 
this context is the realization and control of strongly- 
correlated systems with cold atoms. In this respect, ul- 
tracold dipolar gases play an important role, as the na- 
ture of the dipolar interaction allows one to observe the 
interplay between quantum degeneracy and long-range 
forces. For this reason they are also interesting candi- 
dates for studying statistical theories of quantum long- 
range interacting systems [3|- Dipolar gases of ultracold 
atoms have been experimentally realized with Chromium 
atoms [3, H| • In these experiments the strength of the s- 
wave scattering interaction is conveniently controlled by 
tuning the Feshbach resonance and can be made van- 
ishing, thus leading to the realization of purely dipolar 
systems. Stability against collapse is warranted by po- 
larizing the gas in a two-dimensional geometry, such that 
the dipolar interactions are purely repulsive [3, 0] • 

Different phases of dipolar gases of atoms or polar 
molecules have been theoretically predicted as a func- 
tion of density and dimensionality. In two dimensions, in 
presence of periodic potentials the effect of long-range 
forces gives rise to the appearance of novel quantum 
phases 0|. In a bulk system, recent theoretical work 
predicted the creation of self-organized structures in two 
dimensions: the ground state may exhibit the typical 
features of a crystal or of a quantum fluid depending on 
the density 0, 0] • In a one-dimensional geometry, Lut- 
tinger liquid models describe the long-range properties of 
dipolar gases at ultralow temperatures [9|. In this case, 
by tuning the density, the phase of the gas undergoes a 
crossover from a Tonks-Girardeau gas to a crystal-like 



phase [10[ . The effect of transverse quantum correlations 
in two-dimensional and one-dimensional systems, arising 
from the long-range dipolar force, have been analyzed in 
stacks of pancake traps [~Ll| and in planar arrays of one- 
dimensional tubes [12 ]. Quasi-ordered systems of polar 
molecules in one- and two-dimensions have been recently 
considered for quantum information processing pjj Q • 

One dimensional systems can be exp erimentally real- 
ized with highly anisotropic traps |l5l fl6| . Classically, 
trapped ultracold dipolar systems may exhibit long-range 
order in one dimension at vanishing temperatures, such 
that the equilibrium configuration is solely determined 
by the interplay between the classical repulsive potential 
and the external confinement. In the quantum regime the 
effect of zero point motion is expected to modify sub- 
stantially the crystalline properties. This situation must 
be compared with trapped ultracold ions, interacting via 
Coulomb repulsion, where structural phase transitions 
in selforganized structures have been extensively stud- 
ied pi|[l3]- H ere ? the effects of quantum degeneracy [l8[ 
are negligible in typical experimental setups [l9[ , and the 
structures at experimentally accessible low temperatures 
are essentially determined by the classical potential. 

In this article we study theoretically the competition of 
long-range interactions and quantum fluctuations in the 
transverse stability of a one-dimensional ultracold dipo- 
lar gas in a highly anisotropic trap. We determine the 
properties of the quantum ground state using quantum 
Monte Carlo methods. In particular, we analyze trans- 
verse correlations in the parameter regime in which the 
gas is in the quasi-ordered phase, and where the one- 
dimensional structure is unstable with respect to increas- 
ing the atomic density and/or to reducing the transverse 
potential. We see that transverse patterns are formed, 
giving rise to mesoscopic quantum structures, where the 
transverse density distribution exhibits first two peaks, 
and then, by further opening the trap, multiple peaks. 
Differing from ionic Coulomb crystals, here quantum fluc- 
tuations smoothen the transition from the single- to the 
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double-peaked distribution. Nevertheless, we observe a 
relatively sharp transition to the various transverse struc- 
tures for parameter regimes which can be identified by 
means of a classical theory. 

This work is organized as follows. In Sec. [II] we in- 
troduce the theoretical model for a two-dimensional gas 
of dipolar bosons in presence of tight transverse confine- 
ment and identify the relevant length scales. In Sec. IIIII 
we discuss the phase diagram of the system, which we ob- 
tain using a full quantum numerical simulations; we iden- 
tify the parameter regime where the linear-zigzag chain 
transition could be observed. In Sec. IIVI the conclusions 
are drawn. 



II. THEORETICAL MODEL 

We consider a system of ultracold bosons (atoms or 
polar molecules) of mass m possessing large dipole mo- 
ments, which are confined in the x — y plane, while their 
dipolar moments are aligned perpendicularly to the plane 
by an external field. In this limit the interaction is the 
repulsive dipolar potential 



V(p) 



Cdd 
4ir\p\ : 



where p is the polar coordinate on the plane and Cdd 
denotes the dipolar interaction strength. Here, s-wave 
scattering is neglected, assuming that it vanishes due to 
a properly tuned Feshbach resonance, see for instance 
Ref. 51. The gas is assumed to be homogeneous along the 
x-direction, where it is characterized by the linear density 
n, and it is confined along the ^/-direction by a harmonic 
potential of frequency v t . This external potential sets the 
characteristic length 



O'ho 



h 



determining the transverse size of the single-particle 
wavepacket, and which we choose as unit length. De- 
noting by Qj = pj I dho — (xj , jjj ) the rescaled coor- 
dinates, the corresponding dimensionless Hamiltonian 
H = H/ hvt reads 



~ 1^1 d 2 d 2 „ 2 ^ 



^0 



,|3 



(1) 



where the parameter fo = r^/aho is the characteristic 
length of quantum coherence in dipolar gases, 



TO 



mCdd 
4irh 2 ' 



in units of a^o- 

The effective one-dimensional Hamiltonian is recov- 
ered from Eq. (pp) when the chemical potential \i is much 



smaller than the level spacing of the transverse oscilla- 
tor, \i <C hv t . In the quasi-ordered phase, for nro ^> 1, 
this corresponds to the inequality E^ D ^ /N <C hvt, where 
Ecl^ = N(nro) 3 ((3)h 2 /(mrg) is the potential energy 
of a classical crystal [lOj, and which leads to the re- 

— 1/3 

lation n <C r In the quantum gas regime, for 

nro 1, the gas is essentially described by a Tonks- 
Girardeau gas [20| and the condition to be fulfilled is 
E^q^/N = 7r 2 h 2 n 2 /(6m) <C hv t , which is equivalent to 

the requirement nctho <C 1. The expression for E^q^ /N 
does not depend on the strength of the dipolar interac- 
tion, as the interparticle distance is large and the poten- 
tial interaction can be neglected. For later convenience, 
we introduce the rescaled density 



n = ncLho, 

such that when nr~o ^> 1 the system is in the quasi- 
ordered phase, while for ft <C 1 it is a Tonks-Girardeau 
gas. 



III. PHASE DIAGRAM 

Using the criteria in Sec. [Ill we can identify different 
phases of the ground state of the dipolar system in a 
phase diagram with axes fo and n, and which is reported 
in Fig. dJ Here, the black dashed line corresponds to the 
curve E( 1D ^ /N = hv t , where the energy per particle of 
the one-dimensional dipolar gas, given in [10], is equal 
to the level spacing of the transverse confinement hu t . 
This line separates the two regimes, where the dynamics 
is essentially one or two dimensional. The short-dashed 
line corresponds to the curve nfo = 1, which separates 
the gas from the quasi-ordered phase. These lines are 
to be intended as indicators, the transition of the gas 
from one phase to another (from one-dimensional to two- 
dimensional; from gas to quasi-ordered phases) being a 
crossover. The behavior deep below the black dashed 
line of this phase diagram has been studied in Ref. @, 
where it was shown that the dipolar system is essentially 
a Luttinger liquid with Luttinger parameter K < 1. Al- 
though there is only one phase as the thermodynamic 
functions are continuous, nevertheless signatures of quasi 
order, due to the dipolar interactions, can be identified 
in the energy and in the structure form factor [Tol ]. 



A. Classical theory at T = 

We now focus on the quasi-ordered regime, where 
nfo > 1, and first study the dipolar structure by discard- 
ing the kinetic term in the Hamiltonian (pQ). We stress 
that this approach gives only a qualitative description, as 
long-range order cannot exist in one-dimension at finite 
temperatures in a dipolar gas, nor at T = in a quantum 
dipolar gas. The predictions of this classical model are 
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FIG. 1: (color online) Phase diagram as a function of the pa- 
rameters fo and h. The dashed black line corresponds to the 
curve E 1D /N = hv t and identifies the two regimes, where the 
dynamics is essentially one or two dimensional. The Tonks- 
Girardeau and the classical crystal limits are explicitly in- 
dicated in the plot (dashed-dotted lines). The short-dashed 
grey line, at nfo — 1, localizes the crossover between gas- 
and solid- like phases. The green solid lines separate different 
phases (linear, zigzag, multiple chain) in the classical system. 
The red circles are results of quantum Monte Carlo calcu- 
lations and correspond to the appearance of a double-peak 
structure in the radial density profile (see text and Fig. |4}. 
The size of the symbols denotes the error bars. 



expected to be valid in the asymptotic limit nfo — > oo 
of a quantum system. Keeping this in mind, we assume 
that the dipoles are located along the x-axis with inter- 
particle distance a = 1/n, such that their dimensionless 
equilibrium positions are Xj = j/h and yj = 0. The 
stability of this configuration requires the linear density 
being smaller than a critical value, n < n c , where 

c ~-l/5 

and hence depends on the trap frequency v t , while £ = 
(8/(93C(5))) 1 / 5 = 0.6078... is a constant. For n > h c the 
equilibrium positions are distributed in different planar 
structures, depending on the value of n. The first struc- 
ture encountered is a zigzag configuration, where now 
yj = (— iyb/2 with b fulfilling the equation 



£ 

k=0 



Fk = 



i 



r(- 



12(nr )n 4 
3/2)<(5 + 2fc)(l 



(2) 



-5-2/e\ 



-h 2k 



r(i + fc)r(— 3/2 - k) 

These equations are found by evaluating the stable equi- 
librium points of the potential in Eq. (pQ) [17]. The ex- 
pression for the displacement from the axis of the chain 
close to the critical value h c reduces to the expression 
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FIG. 2: (color online) Radial width yj (y 2 ) /aho as a function 
of n for the classical dipolar gas at T = and nfo = 3 close 
to linear-zigzag transition. The solid line corresponds to the 
solution of Eq. J2]), which describes the transition from a lin- 
ear chain to a zigzag configuration, the symbols to the results 
of a classical Monte Carlo simulation, the dashed line is an 
expansion close to transition point, Eq. (|3]), the arrow shows 
the transition point. 



which provides the critical exponent of the parameter b 
for the classical phase transition. At larger values of n 
the system exhibits abrupt transitions to more complex 
structures. Similar behaviors are also observed in one- 
dimensional Wigner crystals [H, [2l| . Figure [2] displays 
the transverse width of the chain yj (y 2 ) as a function 
of n obtained with a classical Monte Carlo simulation, 
while the solid line is the solution of Eq. ([2]). The critical 
value at which the classical second-order phase transition 
occurs is indicated by an arrow. The critical parameters 
for the linear phase transition in a classical system are 
shown in Fig. [T] as straight solid lines. 

Thermal and quantum fluctuations modify this behav- 
ior, thereby affecting the form of transverse correlations 
as a function of the linear density or of the transverse 
frequency v t . Their effect is analysed below by means of 
numerical techniques. 



B. Numerical methods 

We resort to Monte Carlo techniques to study nu- 
merically the properties of the dipolar many-body sys- 
tem. The classical system at finite temperatures is 
studied by means of a classical Monte Carlo calcu- 
lation, where we sample the Boltzmann distribution 
Pcls — exp{—E/kBT} by using the Metropolis algo- 
rithm. The quantum ground-state properties are deter- 
mined by means of variational and diffusion Monte Carlo 
methods [22j |. The trial wave function is chosen in a Bijl- 



Jastrow form 



N 



N 



ip T ( Pu ...,p N ) =Y[fi(Pi) n f 2 (\Pj - Pk\), (4) 

i=l j<k 

where the one-body term is fi(p) = exp{— ay 2 /a 2 lo }. If 
the system is in the quasi-one dimensional regime, where 
the transverse oscillator is in the ground state, the vari- 
ational parameter is the same as for the ground state of 
the harmonic oscillator, a = 1/2. Outside the quasi-one 
dimensional regime, the system spreads in the radial di- 
rection and accordingly the value of a is reduced. The 
two-body term is chosen as in Ref. [6[ 



Hp)- 



C 2 exp{-^-^ 
1, 



\p\] 



< \p\ < R par 

Rpar < \P\ < L x /2 

L x /2 < \p\ 



(5) 

where i^o(|p|) is the modified Bessel function of the sec- 
ond kind, and coefficients Ci, C2, C3 are fixed by the con- 
ditions of the continuity of the wave function and its first 
derivative. The parameter R par is free, it varies in the in- 
terval < Rpar < L x /2 and is optimized by a variational 
procedure. When the distance between two particles is 
small, the influence of other particles can be neglected 
and /2O0) is well approximated by the solution of the 
two-body scattering problem (see short distance behav- 
ior of (0)). In this way the divergent behavior of the 
interaction potential does not cause any numerical insta- 
bility. At large distances the two-body term is written 
in symmetric way. This ensures that f f (L x /2) = and 
periodic boundary conditions are properly satisfied. 

The results of classical simulations are obtained using 
an annealing procedure. This means that the classical 
simulation is started at a large temperature, hsT/hut ~ 
1, and then the temperature is gradually reduced to very 
small values, where the classical system is frozen. By 
starting the simulation at a high temperature one is able 
to explore more effectively the phase space and this helps 
to avoid the system finding a local minimum instead of 
a global one. We repeat the annealing procedure several 
times for each set of parameters h and fo thus choosing 
the lowest energy. 



C. Numerical results 

In the quantum system we study the transition from 
a linear to a planar configuration by considering the ra- 
dial spreading (y 2 ) as a function of density and trap fre- 
quency. In particular, when the chain is in the ground 
state of the transverse harmonic oscillator, then (y 2 ) = 
a^ G /2, while it becomes larger as the transverse potential 
is excited. In Fig. [3] the radial spreading for the quantum 
mechanical case is shown as a function of n and compared 
with the classical result at T = 0. Quantum fluctuations 
clearly smoothen the transition, and in particular intro- 
duce a basic uncertainty related to the oscillator length 
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FIG. 3: (color online) Radial width \J (y 2 ) /aho as a function 
of h for nro = 3 in the classical and quantum dipolar gases 
at T = 0. Circles correspond to the results of a classical 
Monte Carlo simulation; blue solid line corresponds to the 
solution of Eq. J2}. For h < 1.2 (i.e. for linear and zigzag 
configurations) both results coincide. A sudden transition to 
another, broader, planar structure is observed at larger values 
of h > 1.2. Squares (connecting line is drawn as a guide to the 
eye) results of quantum Monte Carlo calculations, dashed line 
corresponds to the radial width of the transverse harmonic 
oscillator, \J (y 2 ) = dho/VZ- 



dho- Similarly, the transition is smoothened in the clas- 
sical system at finite temperatures. It is interesting to 
notice that the second value, at which one observes an- 
other discontinuity of the classical radial spreading, is a 
transition to multiple chains and it is expected to be of 
first order [21]. The quantum mechanical behavior at this 
point will be object of future studies. The results of the 
quantum MC simulations corresponds to the circles in 
Fig. [1] and have been determined by the critical value of 
the parameters for which two maxima are formed in the 
radial density profile. We note that the position of the 
transition in a quantum system is close to the predictions 
of a classical theory with the agreement being even bet- 
ter as fo is increased. On the opposite, for small values 
of fo the two-peak structure, characteristic for a zigzag 
chain, disappears due to an increased role of quantum 
fluctuations and the crystal gets completely melted. 

Figure 2] displays the radial density profile at fixed nro 
for various values of n, which are distributed across the 
value h c = 0.408 ... at which the classical phase transi- 
tion linear-zigzag occurs for the chosen value of nro. For 
ft < h c the radial density profile is a single-peaked curve, 
a Gaussian of width a/^, centered at y = 0. This is the 
situation we expect for the quantum linear chain. For 
h > n c , the density distribution becomes double peaked, 
the peaks being symmetrical about y = 0. The density 
in the center is however significantly different from zero: 
this effect originates from zero-point quantum fructua- 
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FIG. 4: (color online) Radial density profile n(y) for fixed 
nfo = 3 and at h — 0.2, 0.4, 0.5, 0.6, 0.8 (decreasing the 
height at y = 0). This choice of the parameters corresponds 
to effectively changing the frequency of the harmonic con- 
finement. The radial confinement is normalized to unity, 
n(y) dy = 1. The solid line corresponds to the single- 
particle wavepacket of the transverse harmonic oscillator and 
is plotted for comparison. For these parameters, the value at 
which the classical phase transition occurs is n c — 0.408. 



tions, which prevent the particles from being localized at 
the minima of the potential energy, which for n > n c is 
a double well potential. At larger values of nro, never- 
theless, one recovers the expected zigzag configuration for 
values of n closer to at h c : in this regime the size of zigzag 
spreading is comparable to the oscillator length. The size 
of the region, where quantum fluctuations are relevant, 
shrinks as the interaction strength, here represented by 
the parameter fo, is increased. These planar structures 
are of mesoscopic nature as the number of chains is small. 
Here, the radial density profile cannot be described by the 
local density approximation, which, instead, is generally 
applicable in a macroscopic system. Figure [5] displays 
characteristic examples of the pair-correlation function, 
showing that for sufficiently large values of the interac- 
tion a zigzag distribution of dipoles is observed. 

Thermal fluctuations in the classical systems give rise 
to a qualitatively similar behavior. The curves here ob- 
tained for the quantum system are qualitatively repro- 
duced by using a classical system at the effective tem- 
perature T e ff = fWt/2KB- 



D. Discussion 

The phase diagram in Fig. [1] can be experimentally ex- 
plored by changing the density n, thereby moving verti- 
cally, and by varying the transverse frequency v t , thereby 
moving parallel to the straight short-dashed line. Typ- 
ical parameters of ID experiments [l5j], with a,ho ~ 35 
nm and 60 — 400 atoms per tube of length 15 — 50/im, 
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FIG. 5: (color online) Two dimensional contour plot of pair 
distribution function (n(0)n(x,y)) in a quantum system for 
fixed value of one-dimensional density nfo = 40. Bright colors 
correspond to high values of the pair distribution function. 
Upper plot, fo = 100 (linear chain); lower plot, fo = 200 
(zigzag chain). 



are f « 0.07, n « 0.04 — 0.3. While 52 Cr has the 
largest value of ro = 2.4 nm among all atom species with 
which condensation has been reached [H , the use of polar 
molecules such as CO, ND3, HCN, CsCl with ro = 5 nm 
— 340 jam would permit one to cover regions of the phase 
diagram up to the classical region (nf ^> 1). 

The formation of the zigzag structure can be exper- 
imentally revealed in the structure form factor. The 
latter, in fact, exhibits an additional peak at the y- 
component of the wave vector k y = 27r/6, where b is 
the distance between the two chains forming the zigzag 
structure. 

An important question is whether the system under 
consideration is superfluid or not. The problem of su- 
perfluidity is very delicate when treating (quasi) one- 
dimensional systems. While different definitions of the 
superfluidity coincide in homogeneous two- and three- di- 
mensional systems, they lead to contradict ive conclusions 
when applied to a one dimensional system. Indeed, the 
excitation spectrum of one-dimensional Luttinger liquids 
touches zero for a non-zero value of momentum, k = 2 7m. 
Thus, such systems are always normal from the point 
of view of Landau criterion. Alternatively, one can use 
the winding number technique [24j to calculate the super- 
fluid fraction. Applied to the exactly known ground state 
wave function (for example Calogero-Sutherland system) 
it predicts a completely superfluid system. In the mean- 
field limit a system with short range interactions is well 
described by the Gross-Pitaevskii equation. It is natural 
to think that in this regime, where the description of a 
pure condensate applies, the system is superfluid. 

A possible way out, proposed in Ref. [25j], is based on 
the calculation of the energy dissipation caused by drag- 
ging a small probe through the system. The conclusion 
is that there is always some dissipation, which ranges 
from negligibly small (in Gross-Pitaevskii limit) to the 
same as in a normal system (in Tonks-Girardeau limit). 
There is a continuous crossover from (quasi)-superfluid 
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to normal system. In our case, the limit of small density 
n — > corresponds to Tonks-Girardeau regime, where 
the system is normal. Indeed, in this regime the wave 
function of the system can be mapped onto a wave func- 
tion of an ideal normal Fermi gas. At larger values of n 
the dipolar system corresponds to a Luttinger liquid with 
even stronger interactions than in the Tonks-Girardeau 
regime. As a result our system remains normal for all 
densities. 



IV. CONCLUSIONS 

To conclude, we have studied the ground state of a 
dipolar gas of bosons at T = in low dimensions starting 
from quasi one-dimensional geometry and opening the ra- 
dial trap, so that a two-dimensional structure develops. 
We do an analytical study of the system in the classical 
limit close to the transition to the zigzag configuration by 
applying Landau theory and compare these predictions 
with numerical simulations using classical Monte Carlo 
methods. We determine the phase diagram numerically 
using quantum Monte Carlo methods, and study in detail 
the effect of quantum fluctuations on the linear- zigzag 



transition. The transition from a one-dimensional to a 
planar configuration occurs with the creation of a meso- 
scopic structure in the transverse direction, which ex- 
hibits the main features of the transition from a single to 
a double and then a multiple chain distribution of dipoles, 
while quasi-order is observed also in the transverse pair 
correlation function at large densities. Such patterns are 
characterized by non-local correlations, which arise from 
zero-point fluctuations, and which may be important re- 
sources for the realization of quantum simulators p], [23[ . 
Moreover, the control of atomic patterns has potential 
applications for nanostructuring processes [26[ . 
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